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Abstract 

Let {P n }n>o be a sequence of monic orthogonal polynomials with re- 
spect to a quasi-definite linear functional u and {Q n }n>o a sequence of 
polynomials defined by 

Qn{x) — P n (x) + S n Pn — 

with t n 7^ for n > 2. 

We obtain a new characterization of the orthogonality of the sequence 
{Qn} n >o with respect to a linear functional v, in terms of the coefficients 
of a quadratic polynomial h such that h(x)v = u. 

We also study some cases in which the parameters s„ and t n can be 
computed more easily, and give several examples. 

Finally, the interpretation of such a perturbation in terms of the Jacobi 
matrices associated with {P n } n >o an d {Qn}n>o is presented. 
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1 Introduction 



Let u be a linear functional denned in the linear space P of polynomials with 
complex coefficients. 

The sequence of complex numbers {it n } rl >o such that u n = (u,x n ), is said 
to be the moment sequence associated with the linear functional u. The corre- 
sponding z-transform X)^o Un z~^ n+1 ^ defines a function S(z) which is analytic 
in a neighbourhood of infinity. This function is known in the literature as Stielt- 
jes function associated with the linear functional u (see [17]). 

The linear functional u defines a bilinear functional on P x P in the following 
way: (p,q) u = {u,pq), for every p, q e P. If H = [(x k , ar)u]fc°;_ denotes the 
Gram matrix of the above bilinear functional with respect to the canonical basis 
of P, then H is a Hankel matrix, i.e. the entries in every antidiagonal are equal. 

The linear functional u is said to be quasi-definite if the leading principal 
submatriccs H n , n > , of H are non singular. In such a situation, there 
exists a sequence of monic polynomials {P n } n >o with degP„ = n such that 
(u,P n P m ) = k n 5 n . m with k n ^ 0. The sequence {P n }n>o is said to be the 
sequence of monic orthogonal polynomials (SMOP) with respect to the linear 
functional u. 

Taking into account that the multiplication operator by a; is a symmetric op- 
erator with respect to the bilinear functional (., .)„, the SMOP {Pn}n>o satisfies 
for n > 0, a three-term recurrence relation 

xP n (x) = P n +i(x) + f3 n P n (x) + , y n P n ^i(x), 

with Po(x) = 1 and P^i(x) — where j n ^ for every n > 1. Conversely, if a 
sequence of monic polynomials {P n } n >a satisfies a recurrence relation like the 
previous one, then there exists a unique quasi-definite linear functional u such 
that {P n }„>o is the corresponding SMOP. This result is known in the literature 
as Favard's theorem (see [?])■ 

For the SMOP {P n }n>o there exists a sequence {v n }n>o of linear functionals 
such that (iVt, Pn) — S nm . This sequence {v n }n>o is said to be the dual sequence 

P u 

of {P n }n>o- Notice that v n = - — , where for every p e P, pu denotes the 

linear functional such that (pu,q) = (u,pq), q G P. The linear functional 
u 

v o = ~, tT is said to be the normalization of u in the sense that (vq, 1) = 1. 

(w,P,f) 

Moreover, the linear functional (x — a)~ 1 u is defined on P by 

Nix, P(x)-P(a), 

(<x - a,y x u,P) = (it, — — P e P. 

x — a 

If T> = (Po, Pi, . . . ) T denotes the column vector associated with an SMOP, 
then the three-term recurrence relation reads xV — J-pV where J-p is a tridi- 
agonal matrix such that the diagonal entries are {[3 n }n>o, the upper diagonal 
entries are 1, and the subdiagonal entries are {7 n }n>i- It is called a monic 
Jacobi matrix (see [7]). 
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Thus, for a quasi-definite linear functional u we have three characteristic 
elements 

(i) The Stieltjes function. 

(ii) The SMOP {P n } n > . 

(hi) The monic Jacobi matrix J-p. 

The connections between these three basic elements in the theory of orthog- 
onal polynomials with Pade rational approximants, quadrature formula, and 
spectral theory of symmetric operators have been very fruitful. 

In a more general framework, in [24] and [25] the author deals with quasi- 
orthogonal polynomial sequences {Q n }n>o of order r with respect to a positive 
Borel measure \i. They are defined by Q n (x) — P n (x) + b n>n -iP n -i{x) + ... + 
b rhn ^ r P n ^ r (x) with b n<n - r ^ 0, where {P n }n>o is the sequence of orthonormal 
polynomials with respect to //. When the polynomials Q n are characteristic 
polynomials of some symmetric tridiagonal matrices with positive subdiagonal 
entries, many properties of their zeros can be deduced from such a formulation. 
Furthermore, a positive quadrature formula with n nodes which is exact for 
polynomials of degree at most In — r — 1, < r < n, is based on the zeros 
of the quasi-orthogonal polynomial Q n and, as a consequence, for a fixed n, 
every positive quadrature formula is a Gaussian quadrature formula for some 
nonnegative measure. 

Many authors have been interested in the study of SMOP's associated with 
some perturbations of linear functionals. In particular, the so called direct prob- 
lems deal with some canonical perturbations of linear functionals. For instance, 
the analysis of the parameters of the three-term recurrence relation when a pos- 
itive definite linear functional is multiplied by a polynomial (Christoffel trans- 
formation) has been done from a numerical point of view in [9] , from the point 
of view of the corresponding monic Jacobi matrices in [5], and taking into ac- 
count the associated Stieltjes functions, in [TH] and [21]. The relation between 
the corresponding SMOP in the case of quasi-definite linear functionals has 
been deduced in [2U] . The addition of Dirac masses and their derivatives to a 
quasi-definite linear functional is called Uvarov transformation (see [14], [19] . 
|21|). They play an important role in the spectral analysis of higher order linear 
differential operators whose eigenfunctions are SMOP (see [12]). They are re- 
lated to bispectral problems and self-similar reductions of differential operators 
([18]). The connection between the corresponding monic Jacobi matrices has 
been studied in [5]. 

On the other hand, given a quasi-definite linear functional u and polynomi- 
als p, q the analysis of the quasi-definiteness of a linear functional v such that 
pv = qu has been done in the framework of the so-called generalized Christof- 
fel formula when u is a positive definite linear functional and q/p is a positive 
rational function in the support of u (see [9]). In such a case the connection 
between the monic Jacobi matrices has been done in [8]. For the quasi-definite 
case, the relation between the corresponding SMOP has been given in [20] . 
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In the general situation, if degp = 1 and q is a constant, necessary and 
sufficient conditions for the quasi-definite character of v have been done in [17] . 
There, the expression of the perturbed SMOP in term of the initial one as well 
as the relation between the parameters of the three-term recurrence relation 
is stated. This perturbation is called Geronimus transformation. The relation 
between the corresponding monic Jacobi matrices is analyzed in [5] and [19] . 
Finally, for the Stieltjes functions associated with u and v the connection is 
given in [2"T] . 

When degp = degq = 1, the relation between the SMOP associated with u 
and v is stated in [2J. 

Eventually, when p(x) = x 2 and q is a constant, in [3] the authors focus the 
attention in the invariancc of the scmiclassical character for such a transforma- 
tion. 

All these cases correspond to the so called direct problems, i.e., how a per- 
turbation in a linear functional concerns to Stieltjes functions, monic Jacobi 
matrices, and orthogonal polynomials, respectively. 

The study of inverse problems, i.e. to find the relation between two linear 
functionals assuming their corresponding SMOP's are related in terms of an 
algebraic expression, that it is said to be of type k-j, Q n (x) + a n ,n~iQn-i(x) + 
k+l{x) = P n {x) + (a;) + ... + (x) 

has been firstly considered by Geronimus in [TT] in the case 1-2. Later on, the 
relation 1-3 has been analyzed in [6] as well as in [13] from the point of view 
of the corresponding parameters of the three-term recurrence relation. In this 
last contribution the authors also deduce the second order linear differential 
equation satisfied by the perturbed sequence of orthogonal polynomials when 
the initial sequence is a classical one. We also point out that in [22] and [23] 
the authors focus the attention in the study of these perturbed families when 
the initial orthogonal polynomial sequences are the Chebyshev polynomials of 
first and second kind, respectively. In particular, the integral representation of 
the corresponding linear functional is obtained. Such a kind of relations appear 
in the framework of coherent pairs of linear functionals related to Sobolev inner 
products. This problem has been pointed out in [15]. Relations of type 2- 
2 have been completely studied in [T]. In a very recent paper [3J, the authors 
solved a problem studied in a particular case by Grinshpun (see |10|). concerning 
the orthogonality of sequences of monic polynomials {Q n }n>o with Q n {x) = 
P n {x) + sP n -i(x) +tP n - 2 {x) where {P n } n >o is a given SMOP. 

In pQ and [16] , the orthogonality of the sequence of monic polynomials Q n 
is given in terms of the some constant sequences. Moreover, these constants 
are those appearing in the polynomial which relates the corresponding linear 
functional associated with {P n } n >o and {Q n }n>o- More precisely, in [16j for a 
case 1-2 the authors proved the following results: 

Let {P n } n >o be an SMOP with respect to a quasi-definite linear functional 
u with recurrence coefficients {f3 n }n>o an d {7n}n>i- We define a sequence 
of monic polynomials {Q n } n >o by Q n — P n + ji n P n -x, n > 1, where fi n are 
complex numbers and fii ^ 0. The sequence {Q n }n>o is an SMOP with respect 
to a quasi-definite linear functional v if and only P*{x\)x-y — (3q + /ii,0) ^ 
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for all n > 1 and /z„ = - ^^It^^% » where a?i - /3 + Mi, 0)}„> 

denotes the co-recursive sequence at level zero corresponding to the sequence 
{P n }n>o (see 0, [H]), and Xi= fa- (i 2 - j^- 

This condition is equivalent to the fact that for n > 1 the parameters /x n do 
not vanish and satisfy 

Pn - A*n+i = xi, n>l. (1.1) 

Moreover, (x — X\)v = (/3 — x-y — H\)u. (See [THl Theorem 2]). This result will 
be used many times in this paper. 

In the same way in [T], the authors studied the relations 2-2 {P n + s n P n ~i = 
Qn + t n Qn-i)- They obtained the characterization of the orthogonality of the 
polynomials Q n in terms of the existence of two constant sequences. Besides, 
these constants a, a are the values appearing in the relation between the two 
linear functionals, i.e. {x — a)u = k(x — a)v, (a, a, k € C). 

A natural question is: Do these constant sequences depend on the coefficients 
of the polynomial relating the two linear functionals for more general cases? 

The aim of our contribution is twofold. First, we find necessary and suffi- 
cient conditions that the sequence {Q n }n>o, with Q n (x) = P n (x) + s n P n _i(x) + 
t n Pn-2{x), must satisfy in order to be orthogonal with respect to a linear 
functional v. Indeed, if u is the linear functional such that {P n }n>o is the 
corresponding SMOP, then there exists a monic quadratic polynomial h(x) = 
x 2 + ax + b such that h{x)v = ku. We prove that in this case, there exist two 
constant sequences and the values of these constants a and b are precisely the 
coefficients and not the zeros of the quadratic polynomial h(x) = x 2 + ax + b 
such that h{x)v — ku. A connection between the corresponding monic Jacobi 
matrices is stated. 

Second, we obtain necessary and sufficient conditions in order to the above 
relation 1-3 can be decomposed in two relations 1-2 and then proceed by itera- 
tion. This is related with a two step Geronimus transformation. In this case, we 
prove that the constant sequences can be decomposed in another two simplest 
constant sequences, associated with the zeros of the polynomial h{x). Also, as 
a consequence, a matrix interpretation using LU and UL factorization is done. 
Finally, when {P n } n >o is a symmetric SMOP, i.e. P n {—x) = {— l) n P n (x) and 
we assume s n = 0, then we prove that the 1-3 relation yields two 1-2 relations. 

The structure of the manuscript is the following. In Section 2 we present 
the basic background and we give a new characterization of the orthogonality 
of the sequence {Q n }n>o (Theorem l2.2l) in terms of some constant relations for 
{/?„}, {7n}, {s n }, and {t n }. They have not yet studied in the literature. In 
Section 3 we deduce necessary and sufficient conditions for the split of a 1-3 
relation in two 1-2 relations. The matrix interpretation of this problem in terms 
of monic Jacobi matrices is done. Some illustrative examples are considered. 
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2 A new characterization of orthogonality 

From now on, {P n }n>0 wm denote an SMOP with respect to a quasi-definite 
linear functional u which satisfies the three-term recurrence relation 

P n+1 (x) = (x - n )P n (x) - 7 n P n _!(z), n > 0, (2.1) 
P (x) = l, P_i(x) = 0, 

where {(3 n }n>o an d {'Jn}n>i are sequences of complex numbers with j n ^ for 
n > 1. 

Given two sequences of complex numbers {s n } n >i and {t n } n > 2 we define for 
n > 1 the sequence of monic polynomials {Q n }n>o such that 

Q n (x) = P n {x) + s n P n ^i(x) + t n P n ^ 2 (x), with t n ^ 0, n > 2. (2.2) 

First, we characterize the orthogonality of the sequence {Q n }n>o- 

Proposition 2.1 Let {P n }n>o o,n SMOP with recurrence coefficients {(3 n }n>o 
and {"f n }n>i- We define a sequence {Q n }n>o of monic polynomials by formula 

{EM, i.e., 

Qn(x) = P n {x) + S n P n -i(x) + t n P n - 2 (x), U > 1, 

where s n and t n are complex numbers with t n ^ 0, for all n > 2. Then {Q n } n >o 
is an SMOP with recurrence coefficients {f3 n }n>o and {7 n }n>i, 7n 7^ 0, if and 
only if the following formulas hold: 

Sn-lln = Sn7«-1 + *n(Al-2 ~ Pn), n > 2 , (2.3) 
i n -l7n = tnln-2, n > 3 , (2.4) 

where the coefficients (3 n and j n are defined by 

Pn = Pn + s n - S„ + l, n > , (2.5) 
7n = Tn + *n - Wl + S n (Pn-l ~ Pn), n > 1, (2.6) 

wiift so = to = ^i — 0. 

Proof. Inserting formula (|2.1[) in (|2.2p and applying (|2.2[) to a;P„(x), we get 

<3„+i(x) = xQ n (x) + (s n+ i - s n - f3 n )P n (x) 
+ (t n+ i -7„ - s n p n ^i - t n )P n -i(x) 

- (s„7„_i + t n (3 n ^ 2 )P n - 2 (x) - t n J n ^ 2 P n ^ 3 (x), Tl > 1 , 

provided we substitute there and xP„_2(x), using again (|2.ip . Now, 

formula (|2.2[) applied to P n (x) and the definition of /3„ (see (12.51) ). yield 
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Qn+i{x) = (x- j3 n )Q n (x) 

+ [(t n +l - In — tn ~ S n (s n +1 - S n — n + /3 n -l)] Pn-lfa) 

- [Snln-1 + t n {s n +l - s n - P n + (3 n -l)\ Pn-2(x) 

- tn'Jn^Pnsix), n > 1 . 

So {Qn}n>o is an SMOP if and only if there exists a sequence of complex 
numbers (7 ra )n>i, with j n ^ for n > 1, such that 

[Ui+l - 7n ~ tn ~ S n (s n+ i - S n — p n + fin-l)] Pn-l(x) 

- [Snln-1 + t n (s n +l ~ S„ - (3„ + Pn-2)] Pn-2(x) - t„7„_ 2 in-3(^) 

= -%Qn-l(x). 

Furthermore, (3 n and 7„ are the coefficients of the three-term recurrence relation 
for Q n . 

Using again ()2.2j) . straightforward calculations yield (|2.3p and (|2.4| . □ 
This result can be also seen in [6]. 
Remark. Observe that the conditions 7™ 7^ 0, Vn > 3 , follow from (|2.4[) . 

If the sequence {Q n }n>o given by (|2.2p is an SMOP with respect to a quasi- 
definite linear functional v, then it is well known that the linear functional u 
is a quadratic polynomial modification of v, that is, h(x)v = ku with h{x) = 
x 2 + ax + b and k a nonzero complex number. 

In the next theorem, we show that the orthogonality of the sequence {Q n }n>o 
can be characterized by the fact that there are two sequences depending on 
the parameters s n ,t ni [3 n ,^ n which remain constant and these constants are 
precisely the coefficients a and b of the polynomial h. 

Theorem 2.2 Let {P n }n>o be an SMOP with recurrence coefficients {(3 n }n>o 
and {"f n }n>i- We define a sequence {Q n }n>o of monic polynomials by the 
formula &2.ty) , i.e., 

Q n (x) = P n (x) + s n P n -i(x) +t n P n - 2 (x), n>l, 

where s n and t n are complex numbers with t n 7^ 0, for all n > 2 and sq = to = 
t\ = 0. Then the following statements are equivalent: 

(i) {Qn}„>o is an SMOP. 

(ii) There exist two constant sequences {A n } n >i and {i?„} n >i such that, for 
n>l, 

A n = [7n+l + tn+l — t n +2 + S n +i(f3 n — (3 n +l ~ S n +1 + Sn+2)] (2-7) 

tn+l 

+ S n+ 1 - Pri-l - fin = a, 
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B n = [jn+1 - t n +2 + S n +x(P n ~ Pn+1 - S n +1 + S„ +2 )] (2.8) 

tn+l 

X [in + t n - t n+ i + s n (s n+ i - /?„)] + t n - 7„_i + (s„+i - f3 n )(s n - n -l) = b, 



where, by convention, we take 70 = 0. Besides, the values 0/71 and 72 in 
have to be different from zero. 

Furthermore, if u and v are the linear Junctionals associated with the se- 
quences {P n }n>o an d {Q n }n>0; respectively, normalized by (u, 1) = (v, 1) = 1, 
then 

(x 2 + ax + b)v = ku. (2.9) 

with k e C \ {0} . 

Proof. For simplicity we write the proof of this theorem in terms of /3 n and 
7 n and then it is enough to use (|2.5p and (|2.6|) . 

(*) (it) : 

If {Qn}n>o is an SMOP, by Proposition O {23]) and (HH hold. Inserting 
in (|2.3p the expression of j n -i in terms of 7 n +i given by (|2.4j) we get for n > 2, 

-7— -7n - Al-2 - /?n-l + S n -2 = — ^Tn+l ~ fln-l ~ fin + s n~l ( 2 -10) 
tn 

and then we get (|2.7p . 

Now, we will deduce (|2.8|) . 

From definition of 7W, using (|2.4[) we have 

7n = T^-7n+2 + tn - £ n +l + S n (f3 n -i - S n -1 + S n - /?„) , 71 > 1 . 
t n +2 



Multiplying in the above expression by ^f n +i/t n+ 



1 

7«7n+l . , S ^S«7n+1 ln+lln+2 . , 7. ,Sn7n+l . / t n \ 

— h(Sn-l-Pn-l)— = T h(s„-/3 n j— h " 1 7„+l . 

tn+1 tn+1 tn+2 Jri+1 \In+l / 

Using (|2.10p . for n + 1 instead of n, in the expression which appears in 

tn+l 

the right hand side of the above formula, for n > 1 we obtain 

7n7n+l . / a \( S nln+1 a 

— h (Sn-1 ~ Pn-l) — Pi 

tn+1 \ Iri+1 

7n+l7n+2 , S v / S„+l7„+2 5 

= 7 h (S n ~ Pn) 7 Pn+1 

t n +l \ t n+ 2 

tn 



tn+l 



lj 7n+l + S n ([3 n -l + S n - S n -1 ~ Pn) ■ 

Besides, from (|2.4|) and according to the definition of 7 n _i we have, for 
n > 1, 

"Tn+l = Jn-1 ~ tn—1 + t n — S„_i(/3„_2 + S n _l — S n -2 — Pn-l) 
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and, therefore, 

7«7«+i 



, i S n 7n+1 
+ (S„_i - Pn-l) 1 — I Pi 



(2.11) 



tn+l \ tn+l 

+ tn-1 - 7n-l - 7n + S„-l(/3„-2 + Sn-1 - s n-2 - Pn-l) 

7n+l7n+2 . , a \ ( s n+l7n+2 S 

= 7 h (Sn ~ Pn) T Pn+1 

tn+l \ 

+ *n - 7n - 7n+l + S„0 n -i + s n - S n -1 - /5 



Thus, by straightforward computations, using (|2.5p and (|2.6I) we obtain 

(**) (*) : 

First, we will prove that (|2.10p and (|2.11[) yield (|2.4[) . Taking into account 
the new expression of s " +1 ^"+ 2 ^ obtained from (|2.10[) written for n + 1 instead 
of n, we can reformulate (|2.11|) 

lnln+1 , , a \( S "7n+1 3 
— h (s„_i - p n _i) — " Pn 

tn+l \ tn+l 

+ tn-1 - 7n-l + S„_i(/3„_ 2 + S n _i - S n _ 2 - /?„-l) 



7n+l7n+2 



n+1 



+ ( S n - f3n 



Sn7n+1 



Pn-l + Sn-1 



n+1 



+ *n - 7n+l + S n ($ n -i + S„ - s„_i - /3„) . 

Using the definition of 7„_i in the left hand side of the above expression and 
simplifying we get 



7n+l 
tn+l 



In + S n (P, 



n °n 



+ Sn-1 — Pn-l) + tn-1 



7n-l + 



7n+l7n+2 
tn+2 



Using again the definition of "f n we have 

7n+l , 



tn+l 



{in + t n ) = J n - 



ln+lln+2 
tn+2 



that is, for n > 1 



7"+l / tn+l ~ 

In - ~ 7"+2 



tn+l 



f-ri+2 



tn+l ~ 

ln-1 + ~ 7n+l ■ 

tn+2 



(2.12) 



Since 70 — ^72 = 0, we deduce j„ — f^in+2, for n > 1, and therefore (|2.4 



holds. Finally, we point out that from (|2.10|) and (|2.4[) we can derive 

To conclude the proof it remains to deduce the relation between the func- 
tionals u and v in terms of the constants a and b. 
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If we expand the linear functional u in the dual basis < - — > of the 

{(v,Qj) 

polynomials {Qj}j>o (see [T!5]) and taking into account (|2.2[) . then 



2\ t 

j>0 



2 

\ - {u,Qj) ( h r , si 

u= l^ l rV2\ Ql v = ~^-^Ql + —Ql 



Introducing the explicit expression of the polynomials Q2 and Q\ given by the 
three-term recurrence relation we obtain 



to 



7i72 



/ 2 , /«172 S ~o \ , 7172 % 5 l72 - . 5 3 

(a; +(— Po - Pi)x + I — po — 71 + P0P1 



u . □ 



Remark. The proof of the theorem has been written in terms of j3 n and 7„ not 
only for simplicity. If we work in a similar way as above but trying to look for 
that all the expressions appearing therein depend only on j3 n and 7„, then we 
can obtain similar formulas for a and b but a problem appears when we want 
to come back. Indeed, instead of (|2.12|) we obtain 



7n-l 



7n , 



ln+l+tn+1— " , «>0, 



and so we can not achieve (|2.4p . The main reason is that the trivial extension 
of (|2.4[) to n — 2 (fi = 70 = 0) allows us to obtain the value of 70 in terms of 72 
as 70 = tx but not conversely, i.e. 72 in terms of 70. 

In the sequel, we present a matrix interpretation of these results in terms of 
the monic Jacobi matrices associated with the SMOP's {P n } n >o and {Q n }n>o, 
respectively. 

Let V = (Po, Pi, . . . ) T and Q = (Q , Qi, . . . ) T be the column vectors asso- 
ciated with these orthogonal families, and J-p and Jg the corresponding monic 
Jacobi matrices. Then, the three-term recurrence relations for such SMOP's 
read xV — J-pV and xQ = Jq Q. 

Next, we will describe a method to find the matrix Jq using the matrix J-p 
and the polynomial h(x) — x 2 + ax + b. 

Taking into account Q = AAV where A4 = (rrikj) is a banded lower triangu- 
lar matrix such that rrik,k = 1 and mjy = for k — j > 2, then xMV — JqMV 
and, as a consequence, JpV = M~ 1 JqMV. Thus, we get 

MJp = JqM. 

On the other hand, from the classical Christoffel formula (see [9]) we can 
express h{x)V using the matrix representation 

h{x)V =MQ 

where TV is a banded upper triangular matrix such that rikk+2 — 1 and n^ j = 
for j - k > 2. Thus h{x)V = J\fMV, and then 

J v +aJp + bI =MM. (2.13) 
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But, from h(x)Q = MAfQ, we get 

Jq + aJ Q + bl = MM. (2.14) 

As a conclusion, we can summarize our process as follows 

Step 1: Given Jp, we find the polynomial matrix Jp + aJp + bl. 
Step 2: From M and ([2~T5]) we find N. 

Step 3: From (|2.14p we obtain the polynomial matrix Jq + aJg + bl. 

Step 4: Taking into account Jq is a tridiagonal matrix, from step 3 we can 

deduce J Ql since (j Q + § 1) 2 = MAT - (b - ^ I. 

A similar analysis for the truncated matrices yields the principal submatrices 
of Jq can be deduced from a rank-one perturbation of the corresponding leading 
principal submatrices of J-p. 

Let denote {V) n = (P ,Pi, P n ) T and (Q)„ = (Q„, Qi, . . . , Q„) T Then 
the corresponding three-term recurrence relations read 

x{V) n = (J v ) n+1 (V) n + P n+1 e n+1 , (2.15) 

X (Q) n = (Ja)n+l(Q)n + Qn+ie„+l , (2-16) 

where the symbol (A) n stands for the truncation of any infinite matrix A at 
level n, that is, (Jp) n +i denotes the leading principal submatrrix of J-p of size 
(n + 1) x (n + 1), and e n+1 = (0, ...,0, 1) T G K n+1 . 
Taking into account 

(Q)n = {M) n+1 (P) n , (2.17) 

and replacing (|2.17|) in (|2.16j) we get 

X {M) n+ i{V) n = (JQ)n+l(M) n +l(P)n + (Pn+l + S n+1 P n + t n +iPn-l)e n + l 
= [(J Q) n+1 {M) n+ l + e n +l( 

s n+l e n +l + tn+1 rf «+l)](^)n + Pn+ie n +l , 

where d n+1 = (0, . . . , 0, 1, 0) T G M™ +1 . 
Thus, from (f27T5|) 

{M) n +l{Jv)n+l — (^s)n+l(A / ()n+i + & n +i (s n+ ie^ +1 + t n+ id^ +1 ) ■ 
So, in order to obtain ( Jq)«+i we proceed in three steps. 
Step 1: Obtain (M) n +i(Jv)n+i- 

Step 2: Subtract to the previous one the rank-one matrix e n +i(s n +ie^ +1 + 
tn+idn+i)- This leads to a matrix G n +i. 
Step 3: (Jq)„+i - Gn+iftM^+i]- 1 . 

In the case of quasi-orthogonal polynomials associated with positive Borel 
measures, this approach has been considered in [24] and [25] for symmetric tridi- 
agonal matrices. Nevertheless, not every quasi-orthogonal polynomial can be 
represented as the characteristic polynomial of a symmetric tridiagonal matrix. 
Indeed, in [25] Theorem 3.3] a sufficient condition for such a representation is 
given. 
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3 Reducible cases 



Let {P n }n>o be an SMOP with respect to a quasi-definite linear functional u, 
and {Q„} n >o a sequence of monic polynomials strictly quasi-orthogonal of order 
two, satisfying the relation (|2.2[) . 

As a consequence of the Proposition 12 . II we have that {Q rl } n >o is an SMOP 
if and only if the parameters s n and t n satisfy, for n > 3, 

tn+l = t n (l - 7——) + In + S n ((3 n -l - n + S n +1 - S„) , 
tn-1 

Sn+1 = -T^7n-1 + Ai ~ /9n-l + «3 + 7^71 - 02 - fil , 
1-n £2 

with some initial conditions. 

In this section we will study two particular situations where the computa- 
tions of these parameters can be simplified. 

3.1 Iterative case 

Notice that if {Q n }n>o satisfies (|2.2p . then the polynomials Q n cannot be repre- 
sented as a linear combination of the two consecutive polynomials P n and P n -i- 
A natural question arises: Can the SMOP {Q n }n>o be generated from {P n } n >o 
in two steps with the help of an intermediate SMOP? In other words: Do exist 
two sequences of nonzero complex numbers {A ra }„>i and {/J, n }n>i such that 
{R n }n>o defined by R„ = P n + fJ, n P n -i is an SMOP and Q n = R n + A„i?„_i? 
The interest of this question is to simplify the computation of the parameters 
s n and t n , in this case via the parameters A„ and \x n . 

If the answer is positive, then we will say that the representation (|2.2p is 
iterative. For these cases we have s n = X n + fi n and t n = A„^ n _i. 

Now, we give a characterization of the iterative representations when both 
{P n }n>o and {Q„}„> are SMOP's. 

Proposition 3.1 Let {P n }n>o and {Q n }n>a be two SMOP's with respect to the 
quasi- definite linear Junctionals u andv, respectively, which are related by 112. 2\) . 
In this case, the linear Junctionals satisfy (x 2 +ax+b)v = {x—xi){x—X2)v = ku. 
We denote by {f3 n } n >o> {7«}n>i and {(3 n } n >o, {7n}«>i the coefficients of the 
three-term recurrence relation o/{P„}„>o and {Q n } n >o> respectively. 

Then, the representation 112.2]) is iterative if and only if there is a sequence 
of complex numbers {/i n }n>i, with fi\ ^ s% and fj, n 7^ 0,n > 1, such that 

C n = f3 n - ^ n+ i = xi (or x 2 ), n > 1 , (3.1) 

and 

Hn+i = s n+ i - , n>l. (3.2) 
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Proof. If the representation is iterative then there exists a sequence of 
complex numbers {n n }n>i, with fi n ^ 0, such that R n = P n + jjL n P n -\ is an 
SMOP. Then, from Theorem 2 in [IB], 

C n = j3 n - Hn+i = Xi (or x 2 ), n>l . 

fJ-n 

Substituting R n in (|2.2[) . we get 

Q n (x) = R„(x) + (s n - fi n )R n -i{x) + [t n - /i„-i(s„ - n n )]P n _ 2 (x), n > 2, 

and, since Q n is a linear combination of R n and R n -i, (|3.2|) follows. 

Conversely, given {/_t n }n>i in the above conditions, from [161 Theorem 2], 
(|3.ip implies that the sequence {R n }n>o defined by R n — P n + fi n P n -i, n > 
1, is an SMOP with respect to a quasi-definite linear functional w such that 
(x — x\)w — k\u where k\ 6 C \ {0}. 

Taking A„ = s n — /i„, n > 1, we have A n ^ 0, n > 1, and t„ = A n /i n _i, ?i > 2. 
So, we can write 

Qn = Pn + (A n + /!„)P n _i + A„/l„_iP„_2 = i?n + A„i?„_i, 71 > 1 . 

This means that (|2.2[) is iterative. □ 

Remark. There is another constant sequence involved in this topic. Indeed, as 
a consequence of the above proof and according to [THl Theorem 1] , we get 

D n = p n - A„ - -r^ii = x 2 , (or xi), n > , 
An+l 

with Ao = 0. Moreover, we point out that, in general, the characterization 
of the orthogonality of a family of polynomials {Q n }n>o defined by (|2.2[) can 
be given (see Theorem (|2.2I0 in terms of two sequences {A n } and {B n } which 
remain constant and these constant values are precisely the coefficients of the 
quadratic polynomial h(x) which appears in the relation h(x)v = ku. But if we 
have an iterative representation those sequences {A n } and {B n } (which in fact 
are very complicated) can be decomposed in two other simplest sequences {C„} 
and {D n } associated with the zeros of the polynomial h(x). More precisely, 
when the formula (|2.2|) is iterative then C n + D n = —A n and C n D n ~ B n . 

Observe that, from (|3.2j) . /i n can be expressed as a finite continued fraction 

tn+l 

= Sn+1 . 

^n 

t 2 

S2 

Ml 

Next, as a sake of example we are going to construct a family of monic poly- 
nomials orthogonal with respect to a quasi-definite linear functional v satisfying 
x' 2 v — x a e~ x and which has an iterative representation. 
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Example. Laguerre polynomials. 

Let {£"} n >o be the sequence of monic Laguerre polynomials orthogonal 
with respect to the positive definite linear functional u defined by the weight 
function x a e~ x with a > 0. We can take the auxiliary polynomials R n = L^ 1 
satisfying R n {x) = L"(x) + (see [7]) and the new polynomials Q n 

such that Q n {x) — R n (x) + A„ i?„_i(x). Then the polynomials Q n satisfy the 
relation (|2.2[) with s n = n + X n and t n = (n — 1) A n , n > 1. It is well known 
that the recurrence coefficients of L"" 1 ,^ > 0, are (3 n = 2n + a, n > and 
7„ = n(n + a — 1), n > 1. 

Using formula for this case, by induction we can derive that, for a > 
and a ^ 1, the values of the parameters A„ in terms of Ai are 

T(a)( a -X 1 ) + (\ 1 -1)^. 
X n = n =r, \' , , (3.3) 

and then v is quasi-definite if and only if 

T(n)T{a)(a - Ai) + (Ai - l)r(n - 1 + a) ^ 0, n > 1. 

Notice that ifaeN\{l} then A„ is a rational function of n, namely, 

_ r(a)(a — Ai) + (Ai — + n — 1) . . . (n + 1) 
r(a)(a-Ai) + (Ai -l)(a + n-2)...n 

If a = 1, then, by induction, we can also obtain, for n > 1 

_ (Ai -!)(! + 1/2+ ■■■ + l/^ + l 
An " " (Ar - 1)(1 + 1/2 + • ■ • + l/(„ - 1)) + 1 ' ^- 4J 

and w is quasi-definite if and only if 

(Ai - 1)(1 + 1/2 + • • • + l/n) + 1 + 0, n > 1. 

The quasi-definite linear functional u is given in terms of the Laguerre linear 
functional by 

v a - Xi _ 1 ! 

— = -pr^-z e x ) + £ 
vq I (a) 

for a > 0. In particular for a > 1, we can write 



«o r(«) 



Next, as in the general case, we also present a matricial interpretation for 
these iterative cases. 

If w denotes the corresponding linear functional for {i? n }„>o, then (x — 
x\)w — k\u where X\ is a zero of h{x) as well as (x — X2)v = where X2 is 
the other zero of h(x). 
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It is well known (see [S]) that 

(x - x x )P n = Rn+x + r n R n , n>0, (3.5) 

as well as 

(x-X 2 )Rn=Qn+l+r' n Q n , U > , 

with r n ,r' n ± 0. 

As before, V = {P Q ,Px, ...) T ,K= (R , Rx,... ) T , and Q = (Q Q , Q u . . . ) T , 
and then the corresponding three-term recurrence relations read as 

xV = J v V, xK = J n K, xQ = J Q Q. (3.6) 

On the other hand, from R n = P n + /i„_P„_ 1 and (|3.5p we have the matrix 
representations 

TZ = LxV, (x - xx)V = UxK (3.7) 

where Lx is a lower bidiagonal matrix with 1 as diagonal entries and Ui is an 
upper bidiagonal matrix with 1 as entries in the upper diagonal. 
Notice that from (|3.6p and (|3.7|) we get 

J v - xj = UxLx (3.8) 

as well as L 1 " 1 (Jtj. — xil)7i — UxR- In other words 

J n -xxI = LxUx. (3.9) 

Proceeding in a similar way 

Q = L 2 K, (x-x 2 )n=U 2 Q 

where L 2 is a lower bidiagonal matrix with 1 as entries in the diagonal and U 2 
is an upper bidiagonal matrix with 1 as entries in the upper diagonal. Thus, we 
get 

Jn~x 2 I = U 2 L 2 (3.10) 

and 

J Q - X2 I = L 2 U 2 (3.11) 

As a consequence we can summarize the process as follows. 

Step 1: Given J-p, from Lx and (|3.8p we get U\. 

Step 2: From we get J n . 

Step 3: Given Jn, from L 2 and (|3.10[) we get U 2 . 

Step 4: From (|3"TT|) we get J Q . 

Notice that this is essentially the iteration of the so-called Geronimus trans- 
formation in the framework of linear spectral transformations of Stieltjes func- 
tions associated with quasi-definite linear functionals (see |26)]. 

Now, as in the general case, we analyze the case of the truncated matrices. 
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The three-term recurrence relation for the sequence {i? n } r i>o reads 

X (TZ) n = (Jn)n+lO^)n + Rn+ldn+l ■ (3-12) 

Since 

{Tl) n = {L 1 ) n+l {V) n , (3.13) 

where (Li) ra +i is a lower bidiagonal matrix, namely, the (n + 1) x (n + 1) leading 
principal submatrix of L%, replacing (|3.13[) in (|3.12p 

x(V) n = [(£l)n+l] _1 [(-/R)n+l(il)n+l + A*n+l e n+l e n+l] {P)n + Pn+l^n+1 ■ 

Thus, from (|2.15|) we get </p)n+l = (^R)n+l(^l)n+l + Mn+l e n+ie^+i , 

i.e. 

("^R.)n+l = (Ll) n +l[(Jv)n+l ~ Mn+l e n+l e n+l] [(-^l)n+l] 1 ■ 

As a consequence, (Jrc)n+i is similar to the matrix 
(Jv) i - Mn+l e w+ie n+ i 

that is, to the matrix obtained when we substract in the entry (n + 1) x 

(n + 1) of the Jacobi matrix ( Jp)„ + i and the other entries remain invariant. 
On the other hand, we have 

(x - xi) (V) n = (C/i)„+i(7e)„ + Rn+ie n +i , (3.14) 

where (t/i)„+i is an upper bidiagonal matrix. 
According to (|2~T5| and (|3~T3l) , (j3~14|) reads 

((J'p) n+1 -XiI rs+ i)(7 ? ) n +P n +ie n+ i = (J7i)„ + i(Li)„ + i(7 , )„+(-P, l+ i+^n+i-Pn)e„+i . 

Equivalently, 

(<^p)n+l ~ Mn+l e n+l e n+l — = (C^l)n+l(-^l)n+l ■ 

But from ([3"T4]) 

[(il)n+l] 1 ((Jll)n+l ~ XiI n+ i)(Li) n+ i — (Ui) n +i(Li) n+ i . 

Thus, 

(Jll)n+1 — Xil n +1 = (Li) n+ i(Ui) n +l ■ 

In a similar way, 

(Jn)n+l — A rl +ie„+ie^ +1 — X2l n +1 = (^2)n+l(i2)n+l ; 

and then 

{Jq)ti+1 — X2l n +1 = (L2)n+l{U2)n+l ■ 
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3.2 The symmetric case 



Assume that the sequence {P n } n >o is orthogonal with respect to a symmet- 
ric linear functional u. Then (3 n = 0,n > 0, and there exist two polynomial 
sequences {Ki}n>0 and {^C}n>o such that for all n, P 2n ( x ) — V n {x 2 ) and 



It is known (see [7]) that {V,i}„>o and {V*} n >o are SMOP's with respect 
to the linear functionals a u and xa u where (a u , x n ) = (u, x 2n ), n > 0. 

If the sequence {Q n }n>o is defined by (|2.2p , in general, the polynomials Q n 
are not symmetric. More precisely, Q n are symmetric polynomials if and only 
if s n = for all n > 1 . 

Again, if {Q n }n>o is an SMOP with respect to a symmetric linear func- 
tional v, then by a symmetrization process there exist {i? n }„>o and {i?*}„>o 
SMOP's with respect to <j v and xo vl respectively, satisfying Q 2n (x) — R n (x 2 ) 
and Q 2n +i{x) = xR„(x 2 ). In this case, from (|2.2[) . we have for > 1 



and the coefficients tin and tm+\ can be computed using the Theorem 2 in [16] . 
Besides, as a consequence of the Theorem 12.21 we get a = and the relation 
between the linear functionals u and v is: (x 2 + b)v — ku. 

Now, we construct some families of monic symmetric orthogonal polynomials 
which have reducible representations but not iterative representations. 

Examples 

Hermite generalized polynomials. 

Let {H£} n >o be the sequence of monic Hermite generalized polynomials 
orthogonal with respect to the positive definite linear functional u defined by 
the weight function \x\ 2fJ - e~ x with fi > —1/2. Consider the sequence of monic 
polynomials {Q n } n >a orthogonal with respect to a symmetric quasi-definite 
functional v satisfying x 2 v — ku. Notice that, now the representation (|2.2|) is 
clearly not iterative because xv is not a quasi-definite linear functional. 

By a symmetrization process, a u is defined by the weight function x^' 1 ^ 2 e~ x 
and we can write H^x) — Lfi 2 (x 2 ), H^n+ii 1 ) = x Ln +1 ^ 2 (x 2 ). 

From (|3.3[) and (|3.4[) we get the explicit expressions of the parameters tm in 
terms of t 2 , for (i > — 1/2, fi ^ 1/2, and fi = 1/2, respectively. 

Otherwise, the relation between the two linear functionals u and v yields 
xa v = ka u . Then the polynomials i?* are orthogonal with respect to x^ 1 / 2 e~ x 
and thus t2n+i = n for all n > 1. Moreover the quasi-definite linear functional 
v is given in terms of the generalized Hermite linear functional u by 



P 2n+1 (x) = xV*{x 2 ) 



Rn{x) = V n (x) + t 2n V n - 1 (x) 



(3.15) 



R* n (x) = v; l (x) + t 2n+1 v:_ 1 (x) 



(3.16) 



v_ _ \i + 1/2 - t 2 

v ~ r( M + i/2) 



x 



\\x\ 2fl e- x ) + 5 
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for n > —1/2. In particular, for [i > 1/2 



V_ _ /i + 1/2 - t 2 2/i _ a 2 t 2 - 1 

rOx + i/2) |x| M-i/2 d °- 



Chebyshev polynomials of the second kind. 

Let {[/ n } n >o be the sequence of monic Chebyshev polynomials of the second 
kind, orthogonal with respect to the positive definite linear functional u defined 
by the weight function (1 — a; 2 ) 1 ^ 2 X(-i.i) [ x )-> with the recurrence coefficients 
Pn = 0, n > and 7„ = |, n > 1. Consider the SMOP {Q n } „>o orthogonal 
with respect to a symmetric quasi-definite linear functional v satisfying (x 2 — 
l)v = ku. 

By the symmetrization process described above, we have {V n } n >o, {Rn}n>o 
are SMOP's orthogonal with respect to the linear functionals a u and a v , re- 
spectively, where a u is defined by the weight function (1 — x) 1 ^ 2 x~ 1 ^ 2 X(o,i){x)- 
Besides, the recurrence coefficients for the polynomials V n are, for n > 1, 
0£" = 72« + 72«+i = 1/2, 7 ,? 1 = 72n 72„-i = 1/16, and /3 CT " = 7l = 1/4. 

Moreover, the relation between the two linear functionals (x 2 — l)v = ku, 
yields (x— l)cr v — ka u . Now, applying in (13.15|) Theorem 2 in [16j . the parame- 
ters tin satisfy (3% u - t 2n+2 - 777 = 1 for all n > 1, i.e. t 2n+2 = ~{\+uk^) an< ^ 
k = -(§ + h)^-. Thus, by induction, we can derive the values of parameters 
t 2n in terms of t 2 ■ Indeed 

_ 1 n(At 2 + 1) - 1 
t2n - 'A (n-l)(4t 2 + 1)-1- (3 ' 17) 

On the other hand, since (a; — l)a;cr^ = kxa u , we can apply in (|3.16[) Theorem 
2 in [15] . Thus, the parameters i 2n +i satisfy /3* au — t 2n+ z — ^ — 1 for all 
n > 1, i.e. f2n+3 = — (5 + TFTr^r) ■ Again, by induction, we can derive the 
values of parameters i 2n +i in terms of £3, and taking into account that t 2 and 
t 3 are related by t 3 = 4737^4^) we have 

, _ 1 2n(4t 2 + 1) + (4t 2 - 1) 

W " "I 2(n-l)(4t 2 + l) + (4t 2 -l)- (3 - 18) 

Moreover, if t% € C\ {0, 1/2, —3/4} then the linear functional w is quasi-definite 
if and only if, for n > 1 

n(4f 2 + 1) - 1 ^ 0, and 2n(4t 2 + 1) + (4i a - 1) ^= Q. 

The expression of v in terms of it is 

- = -(3/4 + t 2 ) (1 - x^-^x^ix) - (*a + V4)<5i - (t 2 + 1/4)<U. 
See also [23]. 
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Chebyshev polynomials of the first kind. 

Let {T n }„>o be the sequence of monic Chebyshev polynomials of the first 
kind, orthogonal with respect to the positive definite linear functional u defined 
by the weight function (1 — # 2 ) _1 ^ 2 X(-i,i) ( x )- The recurrence coefficients are 
n = 0, n > 0, 7„ = j, n > 2, and 71 = 1/2. Consider an SMOP {<5 n }ri>o or- 
thogonal with respect to a symmetric quasi-definite linear functional v satisfying 
(x 2 — l)v = k u. 

Thus, {V n } n >0: {^n}n>o are SMOP's with respect to the linear function- 
als a u and cr„, respectively, where a u is defined by the weight function (1 — 
x)~ 1 / 2 x~ 1 / 2 X(o,i) ( x )- Besides, the recurrence coefficients for the polynomials 
V n are, = 1/2, = 1/16, n > 2, and 7l a " = 1/8. 

Since [x— l)a v = ka u , applying in (|3.15[) Theorem 2 in [TS], the parameters 
tin are 

1 n(2t 2 + 1) - 1 

tin = ; z-, ; .n > Z. 3.19 

4 (»-l)(2«2 + l)-l ' - V ' 

Using the same arguments, from (|3.16p we obtain 

, _ 1 Mfy + 1) + (2< 2 - 1) 

W -"4 2(n-l)(2t 2 + l) + (2t 2 -l)' (3 ' 20) 

where we have used that t^ — j ■ 

Furthermore, if t 2 G C \ {0, 1/2, —1/2} then v is quasi-definite if and only if, 
for n > 1 

n(2t 2 + 1) - 1 ^ 0, and 2n(2< 2 + 1) + (2i 2 - 1) 7^ 0. 
t> is given in terms of u by 

- = ^4^(1 - *V((1 - - 2 )- 1/2 X(-M) (*)) + *i + 
Vq Ztt 

. See also l22l. 
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